SURVIVAL OF ^Li AND ^Li IN METAL-POOR STARS 
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^2 ', Abstract. The relationship between the depletions of ^Li and ^Li is studied 

for two models of lithium burning, below the convective zone. The parame- 
ters of the depletion models are submitted to the constraint that the slope 
O I of the Li theoretical depletion curve agrees with the slope of the observed 

depletion curve, for cool subdwarfs. Other less restrictive models are also 
considered. 

In all cases, a ^Li depletion less than 0.5 dex implies a ^Li depletion less than 
0.1 dex. With the constraint on the slope of the ^Li curve, the depletion of 

K> I ^Lifor the same depletion of ®Liis below 0.05 dex. 

C^ ' The still unsolved problem for the true ^Li abundance in subdwarfs is 

the possible influence of temperature inhomogeneities, raised by Kurucz, 
subsequently shown to be small in the solar case, but not yet computed 
with the inclusion of departure from LTE for metal-poor stars. 



1. Introduction 

Lithium is one of the great fossil records available in our local stellar envi- 
ronment. It was discovered by the Spites (1982) in subdwarfs. 
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Apparently, the ^Li /H ratio, has survived ah the way from 100 seconds 
after the Big Bang, to the birth of the first stars, and, more or less, sub- 
sequently in the atmospheres of these stars, in which we observe it, now, 
13 Gyr later. The value of this ratio, based on the analysis of 40 metal- 
poor stars (Bonifacio & Molaro (1997)), is about 1.6x10"^'^ (by number, 
or 1.1 X 10"^ by mass), close to the Standard Big Bang Nucleosynthesis 
prediction (Schramm & Turner (1998)), for the rj ratio (baryons/photons) 
set by the other products of the SBBN. 

The debated question is that of the survival of ^Li in the atmospheres of 
the stars, during the 12 or 13 Gyr they have spent on the main sequence. 
Two processes can alter the initial lithium content: gravitional settlling at 
the bottom of the convective zone (CZ hereafter) (Vauclair & Charbon- 
nel (1998), and poster at this meeting), and mixing of the matter in the 
CZ with deeper layers , in which the temperature allows nuclear burning 
of Li, (Schatzman (1977), Pinsonneault et al. (1992), Charbonnel et al. 
(1992), Montalban & Schatzman (1996)). Furthermore, Kurucz (1995) has 
suggested that the abundance of lithium determined with plane-parallel 
atmospheres might be strongly in error, due to effects of temperature in- 
homogeneities in the convective zone, closer to the surface in metal-poor 
stars than in solar composition stars. 

We shall discuss here only the main process, the burning process by 
mixing of the CZ with deeper layers, clearly acting on metal-poor stars 
of effective temperature below 5700 K (see fig. 1), and responsible of a 
depletion of ^Li by a factor of 200 in the Sun. 

As ®Li has been also observed in two turnoff metal-poor stars, HD 84937 
and HD 338529, and burns at a lower temperature than ^Li , it is interest- 
ing to see what constraints the observation of ^Li brings on the maximum 
depletion of ^Li . 

2. Burning of ^Liand ^Li during pre-main-sequence 

The burning of these elements during pre-main-sequence has been com- 
puted using the CESAM code (Morel (1997). Opacities are from Rogers 
& Iglesias (1992), Rogers et al. (1996), Alexander & Fergusson (1994), or 
Kurucz (1995) The EFF equation of state is used (Eggelton et al. (1973)). 
The results are given in columns 3 and 6 of table 1. 

It must be noted that there is in all cases some depletion of ®Li during 
the PMS, and that this depletion reaches 0.4 dex at metallicity -1.5 for 
a mass of 0.80 Mq . PMS depletion of '^Liis negligible in the range of 
masses considered here. Depletions are counted as positive for a reduction 
in lithium abundance (e.g. a depletion of 1.0 dex means that the present 
abundance is one tenth of the initial abundance) . 
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TABLE L Depletion of Li and Li, in dex, during the PMS and MS phases ( values 
between parentheses are extrapolated and less reliable) 
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3. Burning of ^Liand ^Li during MS lifetime and after 

3.1. PREDICTIONS BASED ON A SIMPLE CIRCULATION MODEL 

We shall use first a very simple model of destruction based on the fact 
that the rate of destruction of ^Liand '''Li varies as a function of depth 
almost as a step function. In the convective zone itself, and just below, the 
destruction rate is so small that no depletion of these elements occur in the 
time available (12 to 14 Gyr). On a short length, of about one hundredth 
of the radius of the star, the rate increases to such a large value that the 
element is burnt everywhere below, in a time much smaller than the age of 
the object. The depletion in the well-mixed convective zone depends only 
upon its exchange of matter with the zone in which the element is totally 
burnt. The corresponding level is of course deeper for '''Li than for ^Li . For 
example, for a star of mass 0.80 Mq and [Fe/H] =-2.0, the transition occurs 
for ^Liat a depth of 0.155 below the bottom of the CZ, in units of the 
stellar radius, but only at a depth of 0.10 for ^Li . Under the assumption 
that the exchange of matter is driven by convective plumes penetrating 
the radiative zone, as suggested by hydro dynamical simulations (Stein &; 
Nordlund (1989), Rieutord & Zahn (1995)) , the critical parameter for the 
exchange of matter between the CZ and the burning zone is of course the 
separation between the two levels. 

On fig. 1 we have plotted depletions of ^Li according to two recents pa- 
pers (Bonifacio & Molaro (1997), Ryan & Delyiannis (1998)), as a function 
of log(Tejj). We have also computed the distance d of the top of the burn- 
ing zone to the bottom of the convective zone, on the main sequence, with 
the CESAM code, for a grid of models , with masses 0.70[0.05]0.85 Mq ( 
the value between brackets is the step), metallicities [Fe/H]= -2.0[0.5]-1.0 
an helium abundance near Yp = 0.23, and an enhancement of a-elements 
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Figure 1. Abundance of Li as a function of logTe//. The sources are Bonifacio & 
Molaro (1997) and Ryan & Deliyannis (1998). The abundances are given in the usual 
logarithmic scale, in which the abundance of hydrogen is 12.0. These abundances are the 
total abundance of Li, in practice, they are abundances of ^Li , ® Li representing at the 
most 5 per cent of the abundance of '^Li , and most of the time much less. 



by 0.4 dex. As Tgjj is a result of the computation, this ahows to redraw the 
depletion of ^Li as a function of d. In this transformation we have allowed a 
shift of 0.01 between the theoretical log(Te//) and the observed log(re//), 
as shown necessary in Cayrel et al. (1997). 

The result, shown in fig. 2, confirms the idea that d is a key parameter 
for the control of depletion by burning. The two straight lines for the two 
metallicities [Fe/H] = -1.5 and -2.0 can be approximated by the single 
relation: 

D{^U) = 0., for d > O.IGR^ (1) 

i^('Li) 



-2.7{log{d/R^) + 0.83) for d < O.WR^ 



(2) 



The physical interpretation of this is that if the rate of exchange of 
matter between the convective zone and the burning zone is otlj, in fraction 
of mass of the CZ per Gyr, the depletion factor of ^Li after a time t will be: 



R = exp(— aTLjt) 

as the matter returning from the burning zone has a zero-content in Li . In 
decimal exponent the depletion becomes D = 0.4343 *a7Ljt or, for an age of 
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Figure 2. Depletion of ^Li as a function of the distance d between the bottom of the 
convective zone and the top of the burning region. The unit of length is the radius of the 
star. The symbols are : squares: metallicity [Fe/H]=-2.0; circles: [Fe/H]= -1.5; triangles: 
[Fe/H]= -1.0. The straight line is a global fit for the two lowest metallicities. 

13 Gyr: D = 5.65 x otl;. For a depletion of 0.5 dex this gives otl; = 0.088, 
or a circulation rate a little below one tenth of the mass of the CZ per Gyr. 
Because the circulation rate as a function of depth is independent of the 
element considered, formulae (1) and (2) are applicable as well to ^Li, at 
the only condition that d refers to the depth corresponding to the top of the 
burning zone of ^Li and not to the top of the burning zone of ^Li . This is 
how we have computed the predicted values for ^Li , actually derived from 
the observed depletions of Li for the same distance between the bottom of 
the CZ and the top of the corresponding burning zone. The depletions of 
^Liand ^Li obtained are plotted in fig. 3, as open circles. Up to a depletion 
of 1.0 for ^Lithe depletion of ^Li remains below 0.1 dex. The main cause of 
this is the strong depletion of ^Liin the PMS phase. 

3.2. PREDICTIONS BASED ON A RANDOM-MOTION MIXING MODEL 

The first model of mixing was proposed by E. Schatzman (1977), and is gen- 
erally described as mixing by "turbulence", or "diffusion". We shall adopt 
here such a model, phenomenogically described by random-motions char- 
acterized locally by the value of a diffusion coefficient, function of depth. 
If the source of the stirring is convection, one expects the diffusion co- 
efficient to decrease stronly with depth below the bottom of the CZ (see 
Montalban & Schatzman (1996)). Actually we can adjust the value at some 
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Figure 3. Relationship between the respective depletions of ® Li and ^Li . The open 
circles are the results of the simple circulation model of section 3.1. The squares are the 
results obtained with a mixing by random-walk, with a diffusion coefficient adjusted to 
reproduce the depletion curve of ^Li in fig.l. The triangles are for a model with a diffusion 
coefficient independent of depth. 



arbitrary depth and the scale-height of an assumed exponential decline of 
the diffusion coefficient in such a way that the abundance of '^Li follows the 
behaviour seen in fig. 1. In fig. 3 the squares show the relationship between 
the two depletions for this model. The conclusions are the same as for the 
first model. An extremely conservative bound for the depletion of ^Li at a 
given depletion of ^Li is given by a model for which the diffusion coefficient 
is constant with depth (triangles in fig. 3). Even in that case the depletion 
of ^Li is les than 0.1 up to depletions of ^Li of about 0.4. For this last model 
the constraint of reproducing the slope of the depletion curve of fig. 1 is 
not fullfilled. 

4. Conclusions 



Models of destruction by burning , reproducing the slope of the depletion 
curve of ^Li , in unevolved metal-poor stars cooler than 5800 K, show a 
fairly robust ratio of the depletions of ^Liand ^Liin the range of interest. 
The depletion of ^Li remains below 0.1 dex , as soon as the depletion of ^Li is 
below the range 0.5-1.0. This result is valid both for a random motion type 
of mixing or for a stationary circulation type of mixing. A large fraction of 
^Li burning occurs during the pre-main sequence. 

The burning of ^Li leads to an important depletion as soon as one of 
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these two conditions is met: either a mass below 0.80 Mq, or [Fe/H] > —1.5. 
This helps to understand why ^Lihas been found so far only in TO or 
subgiants more metal-poor than -1.5. 

Although there are some indications that other effects, as gravitational 
settling, and temperature inhomogeneities in convective elements are not 
very important for ^Li , a global and consistent treatment including micro- 
scopic diffusion, mixing, mass loss and temperature inhomogeneities is a 
must for the coming years. 
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